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1.  introduction 


Turbulent,  plane  wakes  generated  by  circular  cylinders  in  t lie  absence  of  a  pressure 
gradient  have  been  the  subject  of  numerous  experimental  investigations,  the  most  notable  of 
which  are  those  by  Townsend  (1947,  1949).  The  results  of  these  studies,  along  with  results 
obtained  in  other  free  shear  flows,  led  to  the  early  ideas  of  self-preservation  and  Reynolds 
number  independence  proposed  by  Townsend  (1956).  It  was  postulated  that,  sufficiently  far 
downstream  from  the  cylinder,  an  asymptotic  self-preserving  state  is  achieved  for  which 

the  flow  can  be  described  by  a  single  velocity  scale  and  a  single  length  scale  Lq  (see 
Figure  1.1).  That  is,  the  transverse  distributions  of  mean  velocity  and  Reynolds  stress  must 
be  independent  of  the  streamwise  coordinate  x  when  normalized  by  those  scales.  One 
question  under  investigation  presently  is  to  determine  whether  and  where  these  scales  can 
be  considered  unique. 

It  has  been  well  established  (Townsend  1956,  Uberoi  and  Freymuth  1969,  Symes  and 
Fink  1977,  Yamada  et  al.  1980)  that  far-wake  flows  are  self-preserving  and  that  the 
velocity  and  length  scales,  uQ  and  LQ,  vary  as  (x  -  Xq)'1'/2  and  (x  -  Xq)V2,  respectively,  as 
predicted  by  the  equations  of  motion  and  the  momentum  integral  constraint.  ‘rh.f>  virtual 
origin,  x  j$  used  to  account  for  viscous  or  Reynolds  number  effects.  I;  ;"h  >~nore, 
Townsend  (1956)  indicated  that  at  sufficiently  high  Reynolds  number,  L(y  d  and  'UM  are 
universal  functions  of  x/d  only,  where  d  is  the  diameter  of  the  cylinder.  The  initial 
motivation  for  the  present  study  stemmed  from  comparing  some  early  measurements 
describing  the  downstream  variation  of  the  streamwise  component  of  turbulence  intensity. 
We  observed  large  differences  between  our  data  and  data  reported  in  the  literature  that 
could  not  be  attributed  to  experimental  technique  (Figure  1.2).  Ijt  ;e  dif'eren.:1  s  mn 

the  various  resuits  are  evident,  and  the  trends  in  the  data  are  quite  different  at  large 
x/Cpd,  where  Cq  is  the  drag  coefficient.  At  the  time  icirca  1970),  most  investigators  used 
the  cylinder  diameter  for  the  purpose  of  normalization.  Considerations  bv-d  on  the 
equations  of  motion  show  that  the  momentum  thickness,  9,  should  have  been  used  as  the 


normalizing  length  scale  for  the  small-deficit  wake.  That  is,  the  drag  force  exerted  on  the 
fluid  by  the  cylinder  should  be  used  to  define  the  initial  fiow  conditions.  We  v.iii  use  2'j 
because  Cpd  =  26.  It  can  be  shown  that  the  normalized  velocity  and  length  scales  should 
vary  as 
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The  data  of  Townsend  (1956),  Symes  and  Fink  (1977),  and  Yamada  et  al.  (1980)  confirm  the 
above  relations  and  show  that  mean  velocity  profiles  observed  by  each  investigator  are 
self-similar  when  scaled  by  their  individual  velocity  and  length  scales.  However, 
comparison  of  the  data  indicates  a  possible  lack  of  universality  in  the  behavior  of  these 
scales;  that  is,  different  wakes  developed  at  different  rates  with  downstream  distance. 

Sreenivasan  (1981)  examined  the  manner  in  which  wakes  produced  by  a  variety  of 
generators  approached  self-preserving  states.  He  observed  substantial  differences  in  the 
way  these  flows  evolved,  even  though  each  flow  preserved  the  shape  of  the  mean  velocity 
profile  when  normalized  by  its  own  characteristic  scales.  Sreenivasan  and  Narasimha 
(1982)  suggested  that  a  unique  self -preserving  state  exists  for  all  plane  wakes  and  defined 
the  characteristic  constants  stemming  from  their  suggestion.  We  felt  at  the  time  that  their 
data  did  not  fully  support  their  conclusion,  and  the  present  study  compiles  further  evidence 
negating  it. 

One  possible  explanation  for  the  lack  of  uniqueness  is  suggested  by  the  results  of 
Symes  and  Fink  (1977),  who  investigated  the  effect  of  turbulence  in  the  free  stream  on  the 
development  of  wakes.  They  showed  that  the  relative  scale  of  the  external  turbulence, 
rather  than  the  level  of  the  turbulence  intensity,  was  an  important  parameter  affecting  the 
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development  of  the  wake.  As  each  experimental  facility  has  its  own  unique  free-stream 
disturbances,  wakes  generrted  in  a  particular  facility  may  be  unique  tc  that  ir  nity  orly. 
The  lack  of  uniqueness  of  various  small-deficit  wake  flows  might  be  explained  in  terms  of 
the  instability  of  the  mean  velocity  profile  in  the  wake. 

Flow  visualization  (Cimbala  et  al.  1981)  revealed  that  large  coherent  structures 
develop  far  downstream  from  the  wake  generator  and  these  are  not  necessarily  related  to 
the  vortices  shed  from  the  generator.  This  was  sufficient  evidence  to  suggest  that  the 
wake,  like  the  mixing  layer,  contains  large-scale  coherent  structures  which  may  have  a 
wave-like  behavior.  Traveling,  large  structures  were  observed  in  both  laminar  (Freymuth 
1966)  and  turbulent  (Brown  and  Roshko  1974,  Oster  and  VVygnanski  1982)  mixing  layers,  and 
their  behavior  was  explained  by  an  instability  mechanism  (Michalke  1965,  Caster  et  al. 
1984).  The  similarity  of  the  patterns  occurring  in  both  laminar  and  turbulent  states  is  not 
surprising  in  view  of  the  fact  that  the  instability  mechanism  is  principally  inviscid  and  is 
controlled  by  the  mean  velocity  profiles,  which  are  similar  in  both  situations.  Caster  et  al. 
computed  the  amplitude  distributions  and  the  phase  speeds  of  traveling  waves  associated 
with  large  coherent  structures  in  the  plane  turbulent  mixing  layer,  and  they  obtained  very 
good  agreement  with  experimental  results  by  accounting  for  the  effects  of  mean  flow 
divergence.  The  mean  velocity  profile  in  the  wake  is  also  inviscidly  unstable,  and  its  shape 
is  not  affected  by  transition  from  laminar  to  turbulent  flow,  suggesting  that  a  similar 
analysis  could  predict  the  evolution  of  the  large-scale  structures  in  this  flow  as  well. 

The  stability  and  transition  of  a  plane  wake,  generated  by  a  thin  plate  placed  parallel 
to  a  uniform  flow,  was  considered  by  many  investigators  (e.g.,  Sato  and  Kuriki  1961, 
Mattingly  and  Criminale  1972,  Zabuski  and  Deem  1971).  The  analysis  in  these 
investigations,  however,  was  always  concerned  with  the  immediate  neighborhood  of  the 
trailing  edge,  where  the  velocity  deficit  was  greater  than  60  percent  of  the  free-stream 
velocity.  Sato  and  Kuriki  (1961)  limited  their  analysis  to  a  temporal  evolut  on  of  the 
instability  at  one  location  in  the  flow,  and  Mattingly  and  Criminale  (1972)  considered  the 
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instability  of  the  wake  to  natural  disturbances  in  ix)th  time  and  space  and  concluded  that 
the  spatial  evolution  of  a  traveling  wave  gives  superior  predictions  for  the  disturbance 
characteristics  experimentally  observed.  The  analysis  invariably  assumed  that  the  flow  was 
parallel  and  therefore  was  limited  to  a  prescribed  streamwise  location  in  which  the  width  of 
the  wake  was  defined.  The  parallel  flow  assumption  represents  a  severe  constraint  on 
predicting  growing  disturbances  in  the  wake  because,  in  addition  to  the  local  width  of  the 
flow,  the  characteristic  velocity  scale  must  also  change  as  a  result  of  the  divergence. 

In  contradistinction  to  the  mixing  layer,  the  plane  wake  is  susceptible  to  both 
symmetrical  (varicose)  and  antisymmetrical  (sinuous)  modes  of  instability.  The  varicose 
mode  was  traditionally  disregarded  (e.g.,  Sato  and  Kuriki  1961)  because  calculations  based 
on  the  parallel  flow  approximations  indicated  that  the  most  strongly  amplified  disturbances 
were  sinuous.  Data  obtained  in  this  investigation  attribute  the  lack  of  universality  of  the 
self-preserving  wake,  at  least  partially,  to  the  interaction  between  the  two  modes. 
Certainly,  if  one  is  interested  in  examining  the  wake  over  long  distances,  one  cannot 
neglect  the  varicose  mode  of  instability. 

The  scope  of  the  present  investigation  is  limited  to  fhe  small-deficit  wake  starting 
some  100  momentum  thicknesses  from  the  generator  and  extending  to  1500  momentum 
thicknesses  downstream.  The  flow  was  always  incompressible,  with  frce-stream  velocities 
not  exceeding  35  m/ sec  and  typical  Reynolds  numbers  of  a  few  thousand. 

2.  EXPERIMENTAL  ARRANGEMENT 

The  wakes  were  generated  in  the  University  of  Arizona  wind  tunnel  facility.  The 
30-foot-long  tunnel  is  a  closed-circuit  type  built  by  Kenney  Engineering  of  California  and  is 
nearly  identical  to  its  counterparts  at  Tel  Aviv  University  and  the  University  of  Southern 
California.  The  test  section  is  2  feet  wide,  3  feet  high,  and  20  feet  long.  The  top  and 
bottom  walls,  which  are  adjustable  in  height,  were  adjusted  to  compensate  for  boundary 
layer  growth  and  to  obtain  a  zero  streamwise  pressure  gradient.  A  30-hp,  variable-speed 


motor  with  tachomet  ;r  generator  and  a  motoi  controller  drives  an  axial  flow  *.in  with 
variable-pitch  blades.  The  blades  were  set  to  their  minimum  pitch  angle  to  minim. .:e  large- 
scale  turbulence  generation.  In  this  configuration,  the  fan  easily  supplied  the  necessary 
range  of  speeds  in  the  test  section  required  for  the  present  study,  that  is,  2  to  35  m/sec. 
The  tunnel  is  equipped  with  chilled  water  cooling  coils  just  downstream  of  the  diffuser 
section  and  an  electric  heating  unit  just  upstream  of  the  fan  to  control  the  flow 
temperature.  A  Minco  platinum  resistance  thermometer,  connected  to  a  special  bridge  and 
digital  panel  display  unit,  allows  measurement  of  the  mean  temperature  of  the  flow  to 
±0.05°C. 

The  plenum  chamber  contains  4-inch-thick  hexcell  honeycomb,  five  20-mesh  stainless- 
steel  screens,  and  a  4-foot-long  stilling  section.  Following  this  section  is  the  10:1 
contraction  section  with  a  fifth-order  polynomial  contour  to  insure  separation-free 
acceleration  of  the  flow  to  the  test  section.  The  wake  generators  were  mounted 
horizontally  across  the  2-foot  span  of  the  test  section  at  a  streamwise  location  2  feet 
downstream  from  the  inlet.  Measurements  of  the  velocity  profile  at  this  plane  indicated 
that  the  flow  was  uniform  to  ±0.25  percent.  The  free-stream  disturbance  level  in  the 
streamwise  velocity  component  is  approximately  0.03  percent.  The  free-stream  'peed  was 
monitored  using  a  pitot  tube  placed  1  foot  below  the  wake  generator  and  about  2  inches 
into  the  flow.  The  pitot  tube  was  connected  to  a  MKS  Baratron  pressure  transducer  unit. 

Velocities  were  measured  using  a  rake  of  nine  Disa  55P01  hot-wire  probes  connected 
to  Disa  55M01  and  56C01  constant  temperature  anemometers.  The  rake,  which  was  1.75 
inches  in  total  height,  was  used  to  measure  the  mean  streamwise  component  of  the 
instantaneous  velocity.  The  rake  was  mounted  on  an  internal  traversing  mechanism  with  a 
swept-forward,  thin  extension  arm,  placing  the  probes  upstream  of  any  region  of  flow 
interference  caused  by  the  mechanism.  The  mechanism  permitted  traversing  in  the 
streamwise  and  vertical  directions  with  resolutions  of  0.10  and  0.01  inch,  respectively.  The 
anemometer  signals  were  conditioned  using  buck  and  gain  amplifiers  and  simple  low  pass  RC 
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filters  with  a  6-d3  cutoff  point  at  10  kHz.  The  conditioned  s  ;;mls  .■.or''  s  mt  Oir-'ctly  to 
the  analog-to-digital  converter  in  the  data  acquisition  system.  An  L?i  11/23  data 
acquisition/cn-line  computer  system  was  used  for  digital  signal  processing.  Ti  e  major 
components  of  the  system  include  an  LSI  11/23  microprocessor,  a  15-bit  A/D  converter 
with  10  channels  of  simultaneous  sample-and-hold  circuitry;  a  dual-density,  125-ips,  tape 
drive;  a  160  M-byte  hard  disk;  256  K-byte  static  memory;  a  printer/plotter;  two  Tektronix 
4006  graphic  terminals;  a  Tektronix  4611  hard-copy  unit;  and  a  SKYMNK  array  processor. 
The  multiple  channel  capability  allows  simultaneous,  continuous  sampling  of  up  to  10 
channels,  with  variable  sampling  frequency  up  to  10  kHz.  The  system  can  be  used  as  a  data 
logger,  i.e.t  to  create  digital  tapes,  or  as  an  on-line  computer  for  real-time  analysis. 

For  calibration,  the  hot  wires  were  placed  in  the  free  stream,  well  outside  the  wake, 
along  with  a  pitot  tube  which  was  mounted  near  the  rake.  The  bridge  voltage  signals  and 
the  output  of  the  pressure  transducer  connected  to  the  pitot  tube  were  sent  to  the  A/D 
converter.  An  n^-order  polynomial,  U  =  Pn(E),  where  the  independent  variable  E  is  the 
conditioned  bridge  voltage,  was  fit  to  several  calibration  velocities  for  each  hot  wire, 
thereby  providing  an  overall  calibration  for  each  sensor.  Whenever  the  velocity  deficits 
exceeded  0.20  U  ,  a  feoor.d-order  plyo  vnia!  was  used;  whereas  for  u  <  0.10  U  ,  a  linear 
fit  was  adopted,  speeding  on-line  computations. 

During  the  course  of  the  experiments,  the  temperature  of  the  flow  was  maintained  at 
±0.10°C  of  the  calibration  temperature.  The  hot  wires  were  continuously  coccked  for 
drift.  The  results  were  sensitive  to  any  minor  calibration  changes  because  the  maximum 
velocity  deficit  was  on  the  order  of  5  percent.  Generally,  27  to  36  data  points  were  taken 
to  define  a  mean  velocity  and  turbulence  intensity  profile. 

The  wake  generators  used  in  the  present  experiments  are  described  in  Table  2.1.  The 
circular  cylinders  were  steel,  drill  rods  and  the  screens  and  solid  str  p  were  $t  unless  stesl. 
The  cylinder,  screens,  and  solid  strip  were  all  mounted  under  adequate  tension  to  nsure  that 
the  generators  were  straight  and  rigid.  The  screens  and  solid  strip  were  specially 
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cases.  To  introduce  controlled  excitations  into  the  wake,  a  tma.i  :!.in  .5  n  ’••egtn  was 
attached  to  the  trailing  edge  of  the  flat  plate  and  airfoil.  Violin  strinrs  wo-e  used  to 
connect  the  downstream  edge  of  each  side  of  the  flap  to  two  mat -tv  :  -p  ■■avers,  one 

located  on  each  side  of  the  plate  just  outside  the  tunr.e!  sidewalls.  Tint  speakers  wore 
driven  in  phase  at  the  desired  amplitudes  and  frequencies  by  in  audio  amplifier  fed  by  a 
Krohn  Hite  function  generator.  The  forcing  frequency  and  amplitude  were  monitored  by  a 
frequency  counter  and  rms  meter. 


3.  THEORETICAL  CONSIDERATIONS 

3.1.  Similarity  Conditions 

For  a  developing  wake  flow  sufficiently  far  from  the  generator,  the  transverse 
distributions  of  mean  velocity  and  Reynolds  stresses  are  assumed  to  be  self-preserving. 
That  is,  these  distributions  assume  functional  forms  which  are  independent  of  x  when 
normalized  by  the  velocity  and  length  scjles,  u  end  i()(  respectively.  ibis  can  be 
expressed  in  the  form 


u  =  u.  +  ll0  f(n 
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where  n  -  y/l^  and  u^  anti  L(^  are  defined  n  Figure  1.1.  In  general,  un  and  will  be 


functions  of  the  following  parameters: 
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un,  Ln  =  fcns(x,  P,  u.  U„,  F,  d,  GEO,  uk,.  NFS.  OTHERS. 


'3G  .2) 


where 


=  span  of  the  wake  generator 


=  characteristic  wicfth  of  the  wake  generator 


=  geometry  of  the  wake  generator 


=  amplitude  of  the  free-stream  disturbance  level 


=  nature  of  the  free-stream  disturbance 


=  drag  force  on  the  wake  generator  per  unit  length 


OTHERS  =  magnitude  and  nature  of  any  vibration  of  the  wake  generator 


The  conditions  under  which  self-preserving  flow  is  possible  can  be  obtained  by  substituting 
the  self-preserving  distributions  into  the  equations  of  mean  momentum  and  turbulent  kinetic 
energy  and  examining  the  coefficients  in  the  resulting  equations.  For  the  small-deficit  far 
wake  in  the  absence  of  a  pressure  gradient,  i.e.,  when  Uq/U^,  «  1,  Townsend  (1970) 


obtains  the  conditions 


UcLq  duQ  dt_Q 


dx  uQ  dx 


(3.1.3) 


The  self-preserving  functions  are  also  subject  to  the  momentum  integral 


constraint 


['*.) 


dy  =  a 


(3.1.4) 


where  9  is  the  momentum  thickness.  In  terms  of  the  self-preserving  function,  f,  this 


becomes 
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where 
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f n(  n )  dn 


n  =  1.  2 


(3.1.6) 


J  -o 


are  constants  for  a  given  wake  flow.  For  the  small-deficit  wake,  u_/L'„  «  1,  and 
therefore  the  second  term  in  (3.1.5)  can  he  dropped,  placing  a  constraint  on  the  product 

(uqLq)  of  the  two  scales,  since  in  the  absence  of  a  pressure  gradient,  F,  9,  and  U^,  are 
constants  independent  of  x.  Dimensional  reasoning,  along  with  the  linear  part  of  (3.1.5), 
indicates  the  variables  F,  a,  and  U  in  equation  (3.1.2)  should  appear  in  the  combination 
F/pL)2  (see  also  (3.1.4)).  It  can  be  easily  shown  from  (3.1.3)  and  (3.1.4)  that 


and 


(3.1.7) 
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(3.1.8) 


where  x  the  virtual  origin,  is  commonly  assumed  to  depend  on  the  Reynolds  number  and 
geometry  of  the  generator.  This  indicates  that  0  is  the  appropriate  normalizing  length 
scale. 

If  a  universal  self-preserving  state  xists  independent  of  initial  conditions,  free- 
stream  disturbances,  and  other  parameters,  toe  normalized  oistr  bution  functions  f  and  4'$ 
are  universal  functions  and  u  jnd  |  the  normalized  ity  and  'ngth  '■ca.es,  should 

vary  as 
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where  x  =  (x  -  xq)/2  and  A  and  B  are  universal  constants.  2®  is  used  to  normalize  x  here 
because  29  =  C[)d,  which  is  nearly  equivalent  to  normalizing  by  d  for  the  circular  cylinders 
since  Cq  =  1.  Townsend  (1956)  and  others  used  d  in  presenting  their  cylinder  data. 

According  to  Sreenivasan  and  Narasimha  (1962),  equations  3.1.7  and  3.1.8  may  be 
written  in  the  form 


U°  (1  P/ 

“  U-  l 6  I 


A  =  Iq(x9) 


(3.1.11) 


(3.1.12) 


where  W  and  A  are  universal  constants,  provided  the  small-deficit,  equilibrium  wake  is 
independent  of  init’al  conditions  existing  near  the  generator  (see  also  Narasimha  and  Prabhu 
1972).  These  parameters  are  related  to  the  slopes  in  the  relations  (3.1.9)  and  (3.1.10) 
above.  If  we  define  WQ  and  Aq  by  (3.1.11)  and  (3.1.12),  where  x  is  replaced  by  x  -  xQ, 
then  A  =  2/W5  and  B  =  2A^. 


3.2.  Linear  Stability  Analysis 

The  propogation  of  small-amplitude,  wavy  disturbances  in  a  free  shear  layer  was 
considered  analytically  by  Bouthier  (1972),  Crighton  and  Caster  (1976),  and  Caster  et  al. 
(1984).  Since  the  analysis  applied  to  the  plane  wake  is  identical  to  that  presented  in  the 
latter  reference,  only  the  governing  equations  and  essential  features  will  be  presented  here. 
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The  equation  of  motion  considered  is  inviscid  and  has  the  following  form: 


in 

at 


(3.2.1) 


where  fl  is  the  vorticity  and  U  and  V  represent  the  velocity  components  in  the  x  and  y 
directions,  respectively.  Upon  neglecting  the  nonlinear  terms,  which  are  deemed  to  be 
small,  and  assuming  that  a  given  mean  flow  is  parallel  to  the  first  order  of  approximation, 
the  solution  for  the  perturbation  equation  has  a  general  form: 


<li  =  RP(t>(y)  exp(i(ax  -8t))}  (3.2.2) 

where  RP  stands  for  the  real  part  and  the  eigenfunction  1>(y)  is  defined  by  the  inviscid  form 
of  the  Orr-Sommerfeld  equation: 

[U(y)  -  8/a]($"  -  a2$)  -  u"(y)*  =  0  (3.2.3) 

in  which  the  primes  denote  differentiation  with  respect  to  y.  The  wave  number  n  and  the 
disturbance  frequency  8  are  eigenvalues  determined  by  the  solutions  of  (3.2.3),  which 
decay  exponentially  on  both  sides  of  the  wake.  Namely,  the  boundary  conditions  are 

♦  *(±— )  ±  a4>(±»)  =  0  (3.2.4) 

Although  the  divergence  of  the  mean  flow  may  be  partially  controlled  by  the  stresses 
resulting  from  the  interaction  with  the  disturbances  present,  within  the  realm  of  the  linear 
approximation,  the  slowly  divergent  mean  flow  is  assumed  to  be  prescribed  by  equations 
(3.1.1).  Since  the  mean  flow  is  assumed  to  be  known,  the  conditions  of  self-preservation 
are  not  crucial  to  the  analysis;  the  only  requirement  is  that  the  derivatives  of  the  mean 
stream  function  with  respect  to  x  should  be  much  smaller  than  the  derivatives  with 
respect  to  y  (i.e.,  the  boundary  layer  approximation  applies). 
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B y  analogy  with  the  parallel  flow  problem,  the  perturbation  sclutio..  has  the  form 


rx 


*  =  RP  \A(x)  $(x,  y)  exp  j  «(x)dx  -  &t 


(3.2.5) 


wnere  n(x)  is  a  local  wave  number  and  the  eigenfunction  <j>(x,  y)  changes  only  slowly  with 
x.  Since  one  expects  the  correction  terms  resulting  from  the  slow  divergence  to  be  small 
as  well,  a(x)  and  $(x,  y)  at  a  given  streamwise  location  are  still  determined  locally  by 
equation  (3.2.3),  for  which  the  mean  velocity  field  U(x,  y)  is  known,  and  the  correction 
term  is  defined  by 


A(x)  =  A  exD  .  N(*)  dx 

0  exp  CTT51 


ix 

L  0 


where 
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Mi')  =  { 2a  3  $  ♦  U(  y  -  3a )  -  U”$}*dy 


(3.2.6) 


where  $(x,  y)  is  the  adjoint  function  of  $(x,  y)  given  by 


(U(x,  y)  -  0/al($"  -  a2J]  ♦  2U’$'  =  0 


(3.2.7) 
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with  the  boundary  condition',  presented  by 


q'(±«)  -  5(-J1)  =  0  (3.2.8^ 

When  the  mean  velocity  profile  is  symmetrical  with  respect  to  the  line  y  =  0,  the  inviscid 
Orr-Sommerfeld  equation  admits  both  symmetrical  (varicose)  and  antisymmetrical  (sinuous) 
modes  of  disturbances.  For  parallel  flow,  one  may  take  advantage  of  the  symmetry  and 
substitute  a  boundary  condition  on  the  centerline  for  the  boundary  conditions  on  one  side  of 
the  wake. 


♦(0)  =  5(0)  =  1 

V(  0)  =  5‘(o)  =  o 


(3.2.9) 


for  the  sinuous  mode  of  disturbance,  or 


❖(0)  =  5(0)  =  o 

*'(0)  =  5’(0)  »  i 


(3.2.10) 


representing  the  varicose  mode.  By  virtue  of  the  symmetry,  one  may  usually  confine 
attention  to  the  semi-infinite  interval  (0,  ♦«)  in  the  parallel  flow  computations. 

For  a  given  real  frequency  B  of  the  disturbance,  the  complex  eigenvalues  <*(x)  and 
eigenfunctions  $(x,  y)  and  5(x,  y)  were  evaluated  at  each  of  the  streamwise  locations  of 
interest,  provided  the  imaginary  part  a j( x )  <  0  (i.e.,  the  disturbances  in  the  quasi-parallel 
flow  approximation  are  amplified  in  the  downstream  direction),  and  these  solutions  were 
used  in  determining  A(x).  The  mean  velocity  field  used  in  solving  , equations  3.2.3,  3.2.6, 
and  3.2.7  was  obtained  experimentally  and  was  expressed  by  the  exponential  distribution 
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^-  =  1  -  exp[ -0.637 n 2  -  0.056n4]  ( 1.2.11] 

CO  '  CO 

where  Ua  is  the  free-stream  velocity  and  Uq(x)  and  Lq(x)  are  the  velocity  and  length 
scales  discussed  in  section  3.1.  Despite  the  fact  that  the  normalized  shape  of  the  velocity 
profile,  f(n),  remained  invariant  for  all  wake  generators  considered,  the  eigensolutions  had 

to  be  re-evaluated  for  each  generator  separately  since  uQ(x)  and  LQ(*'  dependent  on 
the  initial  conditions,  even  for  the  small-deficit  wakes. 

Sato  and  Kuriki  (1961)  considered  the  temporal  evolution  of  the  small-amplitude, 
sinuous  disturbance  at  a  single  location  in  the  flow  at  which  (1  -  U^/u^)  -  0.692. 
Mattingly  and  Criminale  (1972)  extended  these  calculations  to  both  (nodes  of  instability 
evolving  spatially  and  temporally  at  five  prescribed  locations  in  the  immediate  vicinity  of 
the  trailing  edge  of  a  flat-plate,  wake  generator  for  which  0.44  >_u  /u„  >0.  Since  these 
computations  are  strongly  dependent  on  u^/u^,  the  solutions  obtained  are  of  little  value  in 
predicting  the  character  of  the  amplified  disturbances  in  the  small-deficit  wake  for  which 
Uq/U,  «  1.  Furthermore,  the  assumption  of  parallel  mean  flow  (i.e.,  the  constancy  of 
uq/U<jo  and  of  L(y/0)  led  to  the  general  belief  that  the  varicose  mode  has  a  neglig.ble  effect 
on  the  flow  because  its  rate  of  amplification  (-a;)  j5  smaller.  It  will  be  shown  later  that 
the  interaction  between  the  two  modes  of  instability  leads  to  physically  acceptable  flow 
patterns  associated  generally  with  large  coherent  structures  contained  in  the  wake,  in  spite 
of  the  presence  of  the  small-scale  turbulence,  which  was  not  considered  in  the  calculations. 

4.  EXPERIMENTAL  RESULTS 

4.1.  Ceneral 

The  mean  flow  field  in  the  self-preserving  reg'on  of  a  wake  and  the  inten.itv  of  the 
longitudinal  component  of  the  turbulent  fluctuations  were  measured  for  a  variety  of  two- 
dimensional  wake  generators.  Data  were  obtained  at  dist  inct':,  ranging  from  I'M)  to  2000 
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mcmentum  thicknesses  downstream  of  the  generator,  whom  tvpical  velocity  deficits  on  the 
centerline  of  the  wake  varied  from  0.15  U  to  0.03  U  .  The  velocity  at  which  the 
measurements  were  done  varied  from  7  to  20  m/sec,  corresponding  to  a  range  of  Reynolds 
numbers  based  on  the  momentum  thickness  and  the  kinematic  viscosity  of  air  of  650  to 
3200.  In  most  instances,  the  shape  and  the  size  of  the  wake  generator  were  tailored  to 
provide  a  constant  momentum  thickness.  The  effects  of  Reynolds  number  and  aspect  ratio 
were  examined  separately  in  wakes  generated  by  circular  cylinders.  Forced  sinuous 
disturbances  were  generated  in  the  wake  of  a  symmetrical  airfoil  (having  a  maximum 
thickness  to  chord  ratio  of  30  percent)  and  a  flat  plate  (thickness  to  chord  ratio  of  2 
percent)  by  an  oscillating,  small  flap  mounted  at  the  trailing  edge.  The  results  are 
summarized  in  Table  4.1. 

4.2.  The  Universality  of  the  Mean  Flow  Field  in  a  Small-Deficit,  Plane  Wake 

The  wake  investigated  most  extensively  in  the  literature  is  generated  by  a  circular 
cylinder,  and  therefore  our  investigation  started  with  this  wake  generator.  Some  400 
momentum  thicknesses  downstream,  the  velocity  scale  uQ  was  indeed  proportional  to  x'1  “ 
and  the  width  of  the  wake  was  proportional  to  x^-,  suggesting  that  similarity  of  the 
mean  flow  was  indeed  attained.  The  mean  velocity  profiles  were  plotted  in  the  similarity 
coordinates  and  collapsed  quite  neatly  onto  a  single  curve  described  by  the  exponential 
function 


f( n )  =  exp(-0.637n^  -  0.056n4]  (4.2.1) 

(Figure  4.2.1).  The  exponential  function  traditionally  used  to  describe  the  mean  velocity 

profile  (i.e.,  f ( n )  =  exp[ -0.693  n-] )  overestimates  the  mean  velocity  measured  at  the  outer 
edges  of  the  wake  and,  therefore,  the  fourth-order  correction  term  W3S  acded.  We 
expected  the  flow  to  be  independent  of  Re  when  all  length  scales  were  normal! zed  by  the 


momentum  thickness,  and  indeed  the  values  of  V*  nrrj  were  not  afreet  >d  by  changes  in 
Re  ranging  from  Red  =  1360  to  Red  “  5<jC0*  rat.o  of  the  c/.indtr  was  varied 

from  96  to  384  by  changing  the  diameter  of  the  cylinder  and  keeping  the  span  constant 
(i.e.,  the  distance  between  the  sidewalls  of  the  wind  tunnel).  The  two-dimensionality  of 
the  mean  flow  field  was  checked  at  x  =  430  by  comparing  velocity  and  turbulence  intensity 
profiles  obtained  at  z/L^  =  0  -  7.5  and  was  found  to  be  satisfactory.  The  values  of  Wq  = 
1.75  and  A^  =  0.287  fitted  quite  well  all  wakes  generated  by  a  circular  cylinder  (the  values 
were  averaged  over  experimental  points  with  x  >  200).  This  result  would  give  credence  to 
the  universal  equilibrium  concept  suggested  by  Sreenivasan  and  Narasimha  (1982),  except 
that  the  asymptotic  values  of  the  constants  suggested  by  these  authors  were  quite 
different,  i.e.,  WQ  -  1.63  and  Aq  =  0.300. 

The  universality  of  the  mean  velocity  profile  generated  by  a  circular  evlinder  in  our 
facility  could  lead  to  the  conclusion  that  either  the  development  of  the  plane,  small-deficit 
wake  is  not  susceptible  to  external  disturbances  or  the  velocity-dependent  external 
disturbances  in  our  wind  tunnel  (like  the  fan  blade  passage  frequency,  vibrations,  etc.)  are 
negligible  within  the  range  of  variables  considered.  In  order  to  check  the  \ audit”  of  either 
statement,  we  proceeded  to  investigate  wakes  produced  a  variety  of  r .vo-carv'nsional 
generators.  The  first  family  of  generators  considered  were  screens  with  solidity  ratios 
ranging  from  30  to  100  percent  (a  thin  metal  strip  placed  normal  to  the  free  stream).  The 
porous  screens  have  numerous  advantages:  (i)  there  is  no  flow  reversal  in  the  vicinity  of 
the  generator  with  all  its  ensuing  experimental  complications;  (ii)  the  porous  screens  do  not 
generate  vortices  in  the  same  way  as  does  the  circular  cylinder  and,  therefore,  their  drag 
should  not  be  as  sensitive  to  Reynolds  number;  (iii)  the  roll-up  of  vortices  in  the  mixing 
layers  generated  in  the  neighborhood  of  the  screens  is,  in  principle,  predictable  and 
dependent  on  the  solidity  (Caster  et  al.  1984);  and  (iv)  the  drag  of  screens  having  a 
different  solidity  can  easily  be  equated  by  altering  the  size  of  the  screens.  It  W3s  decided 
to  maintain  the  momentum  thickness  of  all  four  screens  constant  in  order  to  avoid  any 
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questions  about  the  unknown  effects  of  Reynolds  number  or  aspect  ratio  in  the  development 
of  the  plane,  small-deficit  wake.  The  shape  of  the  normalized  mean  velocity  profile  is 
identical  to  that  obtained  for  the  circular  cylinder  in  all  cases  considered,  yet  the  values  of 
Wq  range  between  1.67  for  the  70  percent  solidity  screen  to  1.88  for  the  limiting  case 
representing  the  solidity  of  100  percent. 

The  value  of  Wq  decreases  with  increasing  solidity  of  the  screens,  provided  the 
porosity  suffices  to  prevent  flow  reversal  in  the  lee  of  the  screen;  however,  the  value  of  Wq 
for  the  solid  strip  (which  is  regarded  as  a  screen  having  100  percent  solidity)  is  higher  than 
the  value  of  Wq  for  the  most  porous  screen  investigated  (solidity  of  30  percent).  This 
result  is  attributed  to  the  observed  alternate  shedding  of  vortices  from  the  two  separation 
points  on  the  strip.  The  determination  of  the  threshold  solidity  beyond  which  an  alternate 
shedding  of  vortices  starts  to  occur  is  not  within  the  scope  of  the  present  investigation. 

In  the  absence  of  periodic  forcing,  the  lowest  Wq  measured  in  this  experiment  was 
attained  in  the  wake  of  a  nonlifting,  thick,  symmetrical  airfoil  section.  The  mean  velocity 
field  in  the  wake  of  the  flat  plate  was  quite  similar  to  the  velocity  field  produced  by  the 
wake  of  a  circular  cylinder  except  for  the  location  of  the  virtual  origin  xQ.  jhe 
introduction  of  periodic  surging  reduced  the  effective  Wq  frcm  1.71  to  1.32  wtnout 
affecting  the  overall  drag  experienced  by  the  body  (i.e.,  the  product  W  &  ).  yl  e  location 
of  the  virtual  origin  moved  downstream  with  increasing  amplitude  of  the  forced  oscillations 
(Table  2.1). 

The  product  WqAq  represents  the  conservation  of  momentum  within  the  context  of  the 
linearization  and,  therefore,  its  constancy  for  all  wake  generators  (W  =  0.507  -  0.002) 
reaffirms  the  suggestion  that  the  normalized  shape  of  the  mean  velocity  profile  (Figure 
4.2.1)  is  universal  for  all  the  wake  generators  considered.  An  example  of  the  variation  of 
(Uw/Uqj2  with  x  is  shown  in  Figure  4.2.2,  which  establishes  that  (U^/u q)2  -  Ax  and  that 
the  coefficient  of  proportionality  A  depends  on  the  wake  generator. 
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The  mean  velocity  in  the  wake  obeys  the  similarity  scaling  tor  x  >  200.  Tu?  different 

slopes  of  the  iir.es  drawn  in  Figure  4.2.2  result  in  a  diversity  of  estimates  for  VQ  (Figure 
4.2.3),  varying  between  1.3  and  1.9  depending  on  the  wake  generator.  The  differences  are 
large  and  are  not  attributable  to  experimental  inaccuracy. 

Only  one  value  of  Was  chosen  for  each  wake  generator,  and  it  had  to  satisfy  the 
dependence  of  both  UQ  and  LQ  on  initial  conditions.  The  choice  of  does  not  affect  the 

dependence  of  or  on  the  initial  conditions,  although  the  absolute  values  of  VV^  and 
may  vary  slightly  if  improper  values  of  x^  were  used.  Figure  4.2.3  demonstrates  this  point. 
The  open  symbols  refer  to  physical  distances  measured  from  the  trailing  edge  of  the 
generator  (i.e.,  xQ  =  0),  and  therefore  the  values  of  WQ  based  on  these  data  define  a 
sloping  curve  rather  than  a  horizontal  line.  This  effect  is  particularly  severe  whenever  x^ 
is  large,  as  it  is  in  the  wake  of  a  circular  cylinder.  It  is  apparent,  therefore,  that  for  the 
range  of  distances  considered,  the  mean  flow  in  the  plane,  small-deficit,  turbulent  wake  is 
dependent  on  the  initial  conditions  set  up  by  the  generator. 

The  preservation  of  momentum  deficit  in  the  wake  produced  by  all  the  generators 
considered  above  is  illustrated  by  plotting  9/L  versus  u„  'U  f  Figure  4.2.4).  For  similar 
velocity  profiles, 


with  In  as  defined  in  equation  3.1.6. 

Equation  4.2.2  describes  a  parabola  (a  solid  curve  in  Figure  4.2.4),  while  the  linear 
approximation  (9/Lq)(U^/Uq)  =  is  drawn  as  a  dotted  line.  The  values  of  7 ]  and  1 2  ere 
2.06  and  1.505,  respectively,  and  are  in  excellent  agreement  with  the  data  of  Sreenivasan 
and  Narasimha  (1982).  Measured  values  of  follow  quite  clearly  the  paralolic  curve 

for  all  Uq/U,,  <  0.14.  The  product  (9/LqKU^/Uq)  is  approximately  1.97  (rather  than  = 
2.06)  and  represents  the  shortcoming  of  the  linearization.  It  is  surprising,  therefore,  that 
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'-'p  and  L0  behave  as  they  do  (u^  ~  and  ~  x* if  spit  t  of  the  tact  that  :..f 

shown  in  Figure  4.2.4  illustrate  a  consistant  deviation  from  hav-.r  appro::. -nose  e. 

4.3.  The  Longitudinal  Component  of  Turbulent  Intensity 

For  a  self-preserving  wake  (Townsend  1956),  the  longitudinal  component  of  the 
turbulent  intensity  u  is  given  by 

G2=U0Sl1<n)  (4.3.1) 

where  g-|i(n)  supposedly  is  a  universal  function  analogous  to  f(n).  Sinc^  the  measured 

distribution  of  g-Ji(n)  consistently  showed  a  larger  degree  of  scatter  than  f(n),  a 

normalized  function  u2/(G2)  —  G(n)  (4.3.2)  was  plotted  in  Figure  4.1.1.  The 

normalized  function  was  used  because  of  the  possibility  that  a  slightly  erroneous  Up. 

resulting  from  drift  in  the  response  of  the  hot  wire,  might  cause  the  apparent  scatter  in 

81l(n)»  The  distributions  of  u^  plotted  in  Figure  4.3.1  were  measured  in  the  wakes 

generated  by  the  airfoil,  the  solid  strip,  and  the  70  percent  solid  screen.  Although  each 

distribution  is  approximately  self-similar  (if  one  is  willing  to  disre.’  lrd  the  exco.oivo  s,_..:tur 

in  the  data),  the  function  C(n)  clearly  depends  on  the  type  of  generator.  A  ge  'd  measure 

in  assessing  the  lack  of  universality  is  the  ratio  u2  /u2  (where  the  subscript  CL  refers  to 

CL  max 

the  center  plane  of  the  wake),  which  varies  from  0.73  for  the  solid  strip  to  0.9  for  the 

airfoil.  This  ratio  is  correlated  with  the  rate  of  spread  of  the  wake  or  the  decay  of  the 

mean  velocity  on  its  centerline.  Indeed,  the  lack  of  universality  of  nre<ented  bv 

Cl 

various  investigators  for  the  circular  cylinder  triggered  the  present  investigation  (Figure 
1.2).  The  peculiar  shape  of  G(n),  which  has  a  local  minimum  at  n  -  u,  is  a  s.  irud  i 
single  row  of  large  vortices  distributed  on  the  centerline  of  the  wake  (see  also  VVygnanski 
et  al.  1979).  It  was,  therefore,  anticipated  that  large  coherent  structures  which  retain 
their  characteristic  shape  and  associated  velocity  perturbation  are  responsible  for  the 
apparent  dependence  of  the  wake  on  the  conditions  at  its  origin. 


•  / 


4.4.  Spectra  and  Normal  Modes  of  Instability 

Spectra  of  the  str:amwise  component  of  the  velocity  fluctuations  were  measured  at 
several  x  locations  in  the  wake  behind  the  flat  plate.  The  measurements  presented  in 
Figure  4.4.1  were  taken  on  the  centerline  at  102  <  x  <  587,  while  the  spectra  in  Figure 
4.4.2  were  measured  along  a  curve  representing  the  outer  boundary  of  the  wake,  i.e.,  at  n 
=  3.  Only  the  most  significant  decade  of  the  spectra  have  been  plotted.  T he  abscissa  on 
these  figures  is  frequency  plotted  on  a  logarithmic  scale,  while  the  ordinate  is  f F ( f )  in  order 
to  represent  the  relative  contribution  to  the  streamwise  component  of  intensity  at  a  given 
frequency  f  [i.e.,  u2  •*  / F(f )df  =  /fF(f  )d(logf )]. 

The  spectra  measured  on  the  centerline  of  the  wake  (Figure  4.4.1)  have  a  shape  which 
is  typically  observed  in  any  unbounded,  turbulent,  shear  flow  f Champagne  1073).  The 
spectral  distribution  is  broad,  and  the  frequency  range  associated  with  the  most  energetic 
eddies  gradually  moves  toward  lower  frequencies  as  the  flow  develops  in  the  downstream 
direction.  The  insert  in  Figure  4.4.1  shows  a  log-log  plot  of  the  normalized  (to  unity) 
spectra  F(f),  The  frequency  has  been  rendered  non-dimensional  using  the  local  length 
scale,  Iq,  and  U^.  The  similarity  of  the  spectral  distribution  indicates  that  the  length 
scales  associated  with  the  energy-containing  eddies  (\  =  U  /f)  are  proportional  to  the 

CO 

width  of  the  wake  lg  and,  therefore,  lg  is  the  appropriate  characteristic  length  scale  for 
use  in  rendering  the  mean  velocity  field  self-similar  for  a  given  generator. 

The  spectral  measurements  at  n  =  3,  shown  in  Figure  4.4.2,  represent  fluctuations 
induced  by  the  passage  of  the  large  turbulent  structures  in  the  wake  since  the  coordinate  n 
=  3  is  located  outside  the  turbulent  interface.  The  insert  in  this  figure  shews  that,  as  with 
the  centerline  spectra,  these  spectra  also  scale  with  the  local  widtn  L  Note  the  sS't  of 
the  spectral  peaks  toward  lower  frequency  with  ncreas.ng  x . 

Equation  (3.2.3)  was  solved  tor  the  prescribed  local  mean  velocitv  (i^ki  and  numerous 
real  frequencies  to  obtain  the  variation  the  spatial  rate  of  amplitude 
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ir, creasing  distance  from  :ne  flat  plat.*.  Ii:**  r  of  to-"  >  compi.*  :t;on;  n:,i  ite  that  the 

maximum  local  arr.oiincatinn  rat?  starts  toward  lower  fro  ncies  with  .ncn-.isr’g  x.  Figure 

4.4.3  shows  the  spatial  amplification  rates  for  different  frequencies.  The  value  of  x  for 
which  a  constant  frequency  line  intercepts  the  x  axis  (i.e.,  <ij  =  g)  corresponds  to  the 
downstream  location  at  which  a  wave  at  that  frequency  has  undergone  a  maximum  of 
amplification  according  to  parallel,  linear,  inviscid  stability  theory.  \  plot  or  these 
intercept  values  versus  frequency,  shown  also  in  Figure  4.4.3,  gives  what  we  will  refer  to 
as  the  expected  predominant  frequency  in  the  wake  versus  distance  from  the  flat  plate. 
The  measured  predominant  frequency  range,  defined  as  those  frequencies  at  n  =  3  whose 
amplitude  is  within  90  percent  of  the  peak  energy  value,  are  plotted  in  Figure  4.4.3  and 
compared  with  the  computed  results.  The  dashed  lines  show  the  measured  predominant 
frequency  range  similarity  determined  for  n  =  0.  The  good  agreement  between  the 
measured  predominant  frequencies  associated  with  the  passage  of  the  large  coherent 
structures  and  the  most-amplified  frequencies  calculated  using  parallel  flow,  linear  stability 
theory  suggests  that  the  large  structures  observed  in  this  flow  may  be  related  to  the  two- 
dimensional  instability  modes. 

4.5.  The  Amplification  of  Imposed  Sinuous  Oscillations 

Two-dimensional  sinusoidal  oscillations  in  the  direction  normal  to  the  mean  motion 
were  generated  in  the  wake  by  the  motion  of  a  small  flap  hinged  to  tl'e  trailing  edge  of  the 
flat  plate.  The  frequency  of  the  imposed  oscillation  was  matched  to  the  expected  most- 
amplified  fluctuations  measured  at  n  =  3  in  the  region  of  interest.  A  typical  power 
spectrum  measured  with  and  without  forcing  is  presented  in  Figure  4.5.1.  Small  -  amplitude 
oscillations  do  not  affect  the  turbulent  intensity  in  the  wake  nor  do  they  affect  the  xh  sp** 
of  the  spectral  distribution.  The  two  spectra  presented  in  Figure  4.'.1  are  almost  identical 
with  the  exception  of  the  peak  corresponding  to  the  frequency  of  forcing. 


In  order  to  be  sure  that  the  instability  mechanism  in  the  wak  »  'c:<  r.,!y  r  so,.rsibie 

for  the  amplification  of  the  imposed  oscillations,  measurement  spectra  were  re-  uc 
taken  at  one  location  in  the  flow  for  a  variety  of  forcing  frequencies  while  all  other 
controlled  variables  were  unchanged.  Since  the  background  turbulence  level  is  not 
affected  by  the  introduction  of  forcing  (Figure  4.5.1),  the  ratio  of  the  amplitudes  between 
the  imposed  wave  and  the  background  exhibits  the  sensitivity  of  the  wake  to  the  imposed 
sinuous  perturbations.  The  normalized  spectrum  in  the  center  of  Figure  4.5.2  shows  the 
ratio  of  amplitudes  for  the  most-amplified  frequency  at  the  location  considered.  When  the 
frequency  of  forcing  was  either  higher  or  lower  than  the  most-amplified  frequency  for  the 
given  location  and  flow  conditions,  the  ratio  between  the  peak  amplitude  and  the 
background  diminished.  In  view  of  the  symmetry  of  the  response  around  the  most-amplified 
frequency  (Figure  4.5.2),  which  was  repeated  at  other  flow  conditions  (and  therefore  other 
frequencies),  the  possibility  of  resonance  of  the  mechanical  flap  system  was  discarded. 

The  amplitude  and  phase  of  an  artificially  excited  sinuous  wave  at  a  frequency 
corresponding  to  f0/U  =  6.4  x  10~3  (f  =  20  Hz)  were  calculated  for  the  wake  of  the  flat 
plate.  The  mean  flow  required  for  the  solution  of  cqs.  3.2.2,  3.2.6,  and  3.2.7  was  provided 
by  eq.  4.2.1  and  Table  4.1  (i.e.,  =  1.57,  =  0.323,  =  40).  These  equations  were 

solved  at  intervals  of  x  =  10  for  a  rectangular  window  s'artlng  at  x  =  250  and  ending  at  x  = 
750  and  at  intervals  of  0.1  between  0  <_y/29  <  40.  For  a  given  frequency  (8  =  2m f  L^(X  = 
750)/Um),  the  eigenvalues  a(x)  and  eigenfunctions  $(x,  y)  and  $(x,  y)  were  evaluated  at 
each  of  the  51  streamwise  locations.  These  functions,  their  derivatives  with  respect  to 
both  x  and  y  coordinates,  and  the  mean  flow  information  were  used  to  determine  the 
correction  term  for  slowly  divergent  flow  A(x)  (eq.  5.2.61.  Only  th'  rfeamwise 
component  of  the  velocity  perturbation  was  measured  and  compared  with  the  calculations. 

A  typical  distribution  of  normalized  amplitudes  of  the  sinuous  mode  across  the  entire 
wake  is  shown  in  Figure  4.5.3.  The  solid  line  represents  computed  values,  while  the 
triangles  represent  experimental  results  at  x  500.  The  data  were  acquired  by  recording 
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thc  velocity,  together  with  the  sinusoidal  signal  activating  the  :lar.  !'•  •  .o'.-my  si-n.it 
was  phase  averaged  over  2LI)  c/cles  of  the  flap  motion,  and  '  e  Fourier  trans  r  n  of  t  ‘  >• 

average  records  provided  the  phase  and  amplitude  estimates  of  the  spectral  elements  of  the 
veolcity  field,  which  were  quite  free  from  the  random,  turbulent  fluctuations  present  in  the 
original  signals.  There  is  little  doubt  that  the  artificially  excited  wave  tram  c  in  be 
described  by  an  inviscid  linear  model,  in  spite  of  the  fact  that  the  flow  s  fuliv  turuuit  nt. 
One  may  note  that  the  measured  amplitude  distribution  in  Figure  4.5.3  is  not  symmetrical 
about  the  centerline;  the  lack  of  symmetry  is  partly  attributed  to  inter:  M>mce  between  the 
probe  holder  and  traversing  mechanism  and  the  large  structures  in  the  wake. 

A  detailed  comparison  between  the  predicted  amplitude  distribution  of  the  forced 
wave  and  measured  amplitude  profiles  at  eight  streamwise  locations  is  shown  m  Figir-* 
4.5.4.  In  the  bottom  part  of  this  figure,  tha  velocities  were  normalize*!  bv  their  respective 
maxima,  which  are  replotted  at  the  top.  The  computed  maximum  amplitude  at  x  =  '00  was 
assigned  the  value  of  unity,  which  is  the  only  floating  constant  in  this  comparison.  The 
predicted  and  measured  lateral  distributions  of  amplitudes  are  in  good  agreement  with  one 
another,  as  is  the  decay  of  the  maxima  wit!:  increasing  distance  from  the  gene*-  »*  r.  It  ..as 
surprising,  at  first,  to  note  that  the  maximum  amplitude  of  the  forced  w..  .•->  actually 
diminished  with  x,  in  spite  of  the  fact  that  the  quasi-parai  el  solutions  baser)  on  eq.  1.2 .2 
would  predict  amplification  (<ij  <  0).  The  reasons  for  the  apparent  anomaly  stems  from  the 
fact  that  aj  is  a  small  negative  number  in  the  range  of  distances  considered,  and  it  n 
outweighed  by  the  shape  of  the  eigenfunction  whose  maximum  amplitude  diminishes  with 
increasing  x  (Figure  4.5.5).  A  similar  observation  for  an  axisymmetric  jet  was  made  bv 
Strange  (1982).  The  value  of  the  integral  of  the  perturbation  amplitude  across  the  wake 
increases  somewhat  with  increasing  x  because  the  widtn  of  the  w.w  i-u-rei-.es.  in  fact, 
the  product  of  ihe  maximum  amplitude  and  the  l*>cal  width  is  nearly  constant  between  x  - 
325  and  x  =  700.  This  result  could  not  have  been  predicted  bv  the  parallel  flow 
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distribution  of  -j2  shown  in  Figure  4.3.1,  although  the  ratio  [(ib)  / r )  1  i-  ^“neraliv 

lower  than  the  corresponding  bread-band  distribution.  Because  the  data  are  no  longer 
phase  locked,  the  actual  amplitudes  plotted  are  much  larger  than  for  the  imposed  wave 
train  (Figure  4.5.3).  Since  the  calculated  amplitudes  do  not  vanish  on  the  centerline  of  the 
wake  as  predicted  by  the  sinuous  mode  of  inviscid  amplification,  another  mechanism  has  to 
be  considered.  A  hint  svas  provided  by  the  fact  that  the  amplification  rates  (-Oj)  predicted 
on  the  basis  of  the  parallel  flow  approximation  were  overwhelmed  by  the  presence  of 
longitudinal  gradients  in  the  mean  flow.  It  seems  plausible  that  the  same  longitudinal  mean 
flow  gradients  may  have  enhanced  the  relative  importance  of  the  varicose  mode,  which 
would  contribute  to  the  amplitude  of  the  fluctuations  on  the  centerline;  if  not,  non-linear 
effects  and  secondary  instabilities  may  have  to  be  considered. 

The  calculating  procedure  outlined  in  section  4.4  was  repeated  for  f9/U  =  1.3  -  10"2 
and  the  appropriate  mean  flow  parameters  (Table  4.1).  The  calculations  were  done  twice; 
once  for  the  sinuous  motion,  then  a  second  time  for  the  varicose  mode.  By  assuming  that, 
to  the  first  order  of  approximation,  one  may  simply  superimpose  the  effects  of  both  modes 
of  instability,  the  resulting  amplitudes  can  be  calculated.  For  the  purpose  of  onmnarison 
between  computations  and  experiment,  it  was  assumed  that  the  initial  amplitudes  resulting 
from  both  modes  are  equaliy  important;  namely,  that  the  maximum  amplitude  of  the 
varicose  mode  is  equal  to  the  maximum  amplitude  of  the  sinuous  mode  at  some  initial  x 
distance  from  the  generator.  The  results  of  these  calculations  are  plotted  in  Figure  4.6.2 
for  370  <_  x  <^640;  the  symbols  in  the  figure  represent  data  calculated  from  experimental 
results  and  filtered  at  f6  / U  =  1.3  -  10"2.  There  is  a  qualitative  agreement  between  the 
predicted  and  measured  'aterai  distributions  of  the  u^  amplitudes,  indicating  that  both 
mod  -s  of  instability  are  probably  present  and  can  give  rise  to  the  peculiar  profile  of  the  u 
fluctuations  the  wake.  The  agreement  between  theory  and  experiment  in  this  case  is  not 
is  good  as  for  the  forced  sinuous  wave,  suggesting  that  either  the  two-dimensional 
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approximation  is  not  of  value  or  the  simple  linear  superposition  of  both  modes  is  ^adequate. 
Another  possible  error  stems  from  the  nonlinear  terms  neglected  in  the  present  context. 

In  order  to  explore  further  the  importance  of  the  interaction  between  the  two  modes 
of  instability,  we  resorted  to  flow  visualization  using  a  smoke  wire.  The  wire  was 
positioned  at  x  =  350,  while  the  camera  was  located  at  x  =  500.  The  smoke  patterns  were 
photographed  at  Ub  =  3.4  m/sec  (ReQ  =  600).  Large  coherent  structures  are  clearly 
visible  (Figure  4.6.3),  even  in  the  absence  of  any  imposed  oscillations.  These  structures  are 
similar  in  appearance  to  the  Karman  vortex  street  because  they  seem  to  be  comprised  of 
vortices  of  alternating  sign  of  vorticity,  which  are  placed  in  a  staggered  manner  on  both 
sides  of  the  wake  centerline.  Therefore,  neither  the  varicose  mode,  which  requires  that  the 
vortices  appear  in  pairs  distributed  symmetrically  about  the  centerline,  nor  the  sinuous 
mode,  which  requires  vortices  whose  center  coincides  with  the  centerline,  dominates  this 
flow.  Although  the  vortices  are  large  and  coherent  and  have  a  prevailing  wavelength  of 
approximately  15  cm,  there  is  sufficient  irregularity  in  their  shape,  size,  and  position  to 
cause  the  spectrum  of  the  induced  fluctuation  at  the  edge  of  the  wake  to  be  fairly  broad 
(Figure  4.4.2).  The  introduction  of  periodic  forcing  did  not  have  a  significant  effect  on  the 
shape  and  scales  of  the  large  eddies  visualized  in  this  manner.  The  degree  of  two- 
dimensionality  of  these  structures  was  first  estimated  by  placing  the  smoke  wire  parallel  to 
the  circular  cylinder  but  displacing  it  from  the  generator  in  the  lateral  direction  in  order 
that  the  smoke  would  not  be  entrained  by  the  wake  before  x  =  300.  The  resulting  photo 
(Figure  4.6.4)  indicates  that  the  large  eddies  have  a  tendency  to  be  two  dimensional, 
although  the  two-dimensionality  is  by  no  means  perfect.  There  appears  to  be  a  variation  of 
amplitudes  along  the  span  of  the  wake,  as  well  as  phase  irregularities. 

Coherence  spectra  calculated  from  u  fluctuations  sensed  by  two  probes  separated  in 
the  spanwise  direction  provide  a  quantitative  measure  of  the  two-dimensionality  of  the 
various  scales.  The  two-point,  cross-correlation  function  for  stationary  random  variables 
Ujf^  t)  and  uj{^+  r_,  t  ♦  t  )  is  defined  as 


ij(X.f  T  )  =  Ui(£,  t)  Uj(x_*  r_,  t  ♦  T  ) 


■4.6.1 > 


and  the  cross-spectrum, 


00 

u  0  =  27  RijOjL*  U  r'>  el2lTfTdt  =  Cjj(r_,  x_;  f)  -  iQij(x_,  r_,  f)  (4.6.2) 


where  Cjj  „  called  the  cospectrum  and  Qjj,  the  quadrature  spectrum.  The  coherence 


spectrum  is  defined  as 


C?j(x,  r;  f)  ♦  Qjj(x_,  r_;  f) 


Cohii-' -  f)  |Jii(x_;  f)  Fjj(i+  l;  o 


where  Fjj(x^ ;  f)  and  Fjj(x_ r_;  f)  are  the  familiar  (one-point)  energy  spectra,  i.e.. 


(4.6.3) 


«?(x)  = 


The  phase  angle  9  j j  can  be  obtained  as 


Fii(x  ;  f)df 


(4.6.4) 


®ij(x_»  L’  f~>  =  tan_1 


l|Qij 


(4.6.5) 


The  coherence  spectrum  is  bounded,  and  its  value  must  be  between  0  and  1.  We  will 
consider  the  component  Coh-j-)(x  ,  AZ;  f),  which  represents  the  degree  of  spatial  correlation 
between  the  Fourier  components  of  u-j(x_,  t)  and  tq  (x  +  z,  t  ♦  t )  at  the  same  frequency, 
where  kjs  a  unit  vector  in  the  z  or  spanwise  direction. 

A  spanwise  rake  of  six  hot-wire  probes  with  A z  spacing  from  1.1  to  2.54  cm  was  used 
for  the  coherence  data.  Measurements  were  taken  in  the  wake  of  the  flat  plate  with  and 
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without  forcing  at  ;;  =  430  and  646,  with  the  rake  located  at  r.  =  0,  0.6,  and  3.  The  lateral 
position  n  =  0.6  corresponds  roughly  to  the  position  of  maximum  u  activitv.  The  wake 
response  .at  x  =  646  and  n  =  0.6  to  various  forcing  levels  is  shown  in  figure  4.6.5  for  two 
A z  spacings,  1.1  and  25.4.  The  values  of  the  coherence  at  the  forcing  frequency  (50.0  Hr) 
and  its  first  and  second  harmonics  as  a  function  of  amplitude  are  presented.  All  amplitudes 
are  scaled  by  the  highest  amplitude  of  forcing.  Where  the  relative  forcing  amplitude  was 
less  than  0.25,  the  coherence  spectrum  responded  in  a  linear  manner.  For  amplitudes  in  the 
range  of  0.25  to  0.5,  the  value  of  the  coherence  at  the  forcing  frequency  is  nearly 
independent  of  Az  in  the  range  investigated,  at  least  up  to  A  z/L^  =  iu.  The  relative 
amount  of  energy  tied  up  with  the  forcing  frequency  can  be  determined  from  the  spectrum 
of  the  velocity  fluctuations.  For  this  case,  the  spectral  peak  at  50  Hz  was  one  order  of 
magnitude  above  the  'background*  turbulent  fluctuation. 

Figure  4.6.6  shows  the  coherence  spectra  measured  at  x  =  646,  n  =  0.6,  and  A  z/l^  = 
0.4,  with  and  without  forcing.  The  forcing  frequency  was  50  Hz,  corresponding  to  the 
expected  predominant  frequency  at  this  x,  and  the  relative  amplitude  of  forcing  was  0.25. 
The  effect  of  forcing  sharply  enhances  the  value  of  the  coherence  at  the  forcing  fr-quency 
to  0.92,  with  little  effect  on  the  rest  of  the  spectrum.  The  correspondin';  data  for 

=  10  are  shown  in  Figure  4.6.7(a),  where  only  the  data  for  the  forced  'mso  are 
presented.  The  coherence  for  the  unforced  case  at  this  separiHon  vanished  at  all 

frequencies.  For  the  forced  case,  the  entire  correlation  is  contained  in  the  spectral  spike 
at  50  Hz,  for  which  the  coherence  is  0.87.  Similar  resuits  were  obtained  at  forcing  levels 
as  low  as  0.05,  where  the  peak  correlation  was  0.21  and  50  Hz.  4t  n  =  3 
[Figure  4.6.7(b)],  a  much  higher  coherence  at  the  forcing  t'requencv  was  measured  (0.98), 
indicating  that  the  large  structure  in  the  wake  must  be  highly  two  dimensional  to  generate 
such  a  result  for  the  large  spanwise  separation,  A  z/L^  =  iq.  por  tho  unforced  c  w,  the 


coherence  is  nearly  zero  from  10  to  260  Hz.  The  results  indicate  that  a  slight  amount  of 
forcing  generates  a  strong  two-dimensional  wave  in  the  wake  at  the  forcing  frequency. 

Is  the  proposed  model  capable  of  explaining  the  iarge  structures  observed  in  Figure 
4.6.3?  For  this  purpose,  some  streak-lines  had  to  be  calculated.  Since  the  emulations 
were  ai  ned  at  a  qualitative  understanding  of  this  phenomenon,  a  parallel  flow 
approximation  was  invoked  locally.  It  was  assumed  that  the  particles  were  uniformly 

released  at  *  =  300,  where  u^/u^  =  0.06.  The  prevailing  wavelength  gave  rise  to  3  =  0.8, 
which  corresponds  to  the  spatially  most-amplified  sinuous  wave  train  in  this  mean  flow. 
Both  u  and  v  components  of  the  perturbation  velocity  were  calculated  across  the  entire 
wake,  and  the  corresponding  particle  paths  were  established  from  the  equations 

J  -  Ulx(t),  y(t),  t);  ^  =  V[x(t),  y(t),  t]  (4.6.6) 

(For  a  detailed  description  of  the  procedure,  see  Michalke  1965.)  The  coefficient  of 

velocity  perturbation  in  all  these  calculations  was  taken  as  0.025. 

Five  of  the  streak-lines  calculated  taking  only  the  sinuous  mode  into  consideration  are 
shown  in  Figure  4.0.0(a).  Tb"  streak-lines  have  a  sinusoudal  pattern  undulating  ah  ut  the 
centerline  of  the  w.-ka.  The  amplitude  of  the  undulations  increases  with  increasing 
distance  from  the  source;  at  large  distances,  most  of  the  particles  congregate  at  the  outer 
edge  of  the  wake. 

The  corresponding,  most  strongly  amplified  varicose  mode  h3s  a  perturbation  frequency 
of  8  =  0.5  and,  therefore,  corresponds  to  5/8  of  the  frequency  of  the  prevailing  sinuous 
mode  ( not  quite  a  subharmonic  frequency).  The  computed  streak- lines  for  ttv*  varicose 

mode  are  shown  in  Figure  4.6.8(b).  The  particles  in  this  figure  congregate  in  lumps,  which 

are  symmetrical  about  the  centerline  of  the  wake.  The  streak-lines  shown  in  Figure 
4.6.8(a)  or  4.6.8(b)  bear  no  resemblance  to  the  observed  coherent  structures  sect  in  Figure 
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Th^  streak-lines  shown  in  Figure  4.6.fl(c'  represent  the  combined  motion  or  both  -vices 
o;  insta!  ility  when  the  coefficient  of  the  velocity  perturbation  is  still  maintained  .it  0.025. 
The  initial  ratio  between  the  amplitudes  of  u  and  v  in  the  varicose  mode  and  the  sinuous 
mode  was  0.7.  (This  number  simply  appeared  because  the  eigenfunctions  were  not 
normalized;  changing  the  initial  ratio  by  a  factor  of  two  made  no  substantial  difference  in 
the  pattern.)  The  initial  phase  relation  between  the  modes  was  assumed  to  be  zero. 
(Changing  this  number  also  had  no  effect  on  the  basic  pattern  at  some  distance  downstream 
from  the  source  location.) 

The  prevailing  wavelength  of  the  streak-lines  [Figure  4.6.8(c)]  still  corresponds  to  the 
prevailing  wavelength  of  the  sinuous  mode,  but  the  presence  of  the  varicose  mode  not  only 
modulates  the  streak-lines  but  also  contributes  to  an  apparent  chaotic  behavior.  Most 
important  is  the  fact  that  these  streak-lines,  when  rcplotted  on  a  larger  scale,  re  wimble  the 
pattern  observed  in  the  smoke  photographs  (Figure  4.6.3).  Namely,  the  particles 
congregate  in  a  staggered  manner  about  the  centerline  and  the  large  eddies  are,  at  times, 
separated  by  deep  incursions  of  ‘potential*  fluid.  Sometimes,  these  incursions  are  narrow 
(marked  by  the  letter  "N*  in  Figure  4. 6.9)  end,  sometimes,  they  are  wide  (marked  by  ‘  .V). 
Therefore,  it  transpires  that  only  the  combination  of  both  modes  can  successfully  describe 
the  flow. 

The  distributions  of  vorticity  for  the  three  cases  considered  in  Figure  4.6.5  were 
calculated  directly  from  the  eigenfunctions  solved,  because  the  vorticity  perturbation 

a»(y)  =  -(U"/(U  -  B/a)]<t>  (4.6.7) 

and  the  total  vorticity 

H(x,  y,  t)  =  -U’  ♦  0.015  RP{u(y)  exp[i(ax  -  St)]}  (4.6.8) 


which  is,  of  course,  periodic  in  time.  The  vorticity  contours  plotted  in  Figure  4.6.10  were 
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calculated  for  t  =  2-r/S  and  for  355  <_  x  <_  335,  as  in  Figure  1,6.9.  (The  shaded  regions 

correspond  to  negative  vorticity.)  By  neglecting  either  the  varicose  mode  [Figure 

4.6.10(a)]  or  the  sinuous  mode  [Figure  4.6.10(b)],  the  resulting  vorticity  contours  appear 
to  be  very  regular.  The  contours  are  either  symmetrical  or  antisymmetrical  about  the 
centerline  and  indicate  intensification  of  vorticity  with  increasing  x.  The  vorticity 
contours  resulting  from  the  combined  two  modes  of  instability  [Figure  4.6.10(c)]  are 
surprisingly  similar  to  the  contours  produced  by  the  sinuous  mode  alone.  One  may  notice, 
however,  that  the  intensity  of  the  contours  in  Figure  4.6.10(c)  are  modulated  by  the 

varicose  mode,  even  if  they  are  not  severely  distorted  by  the  addition  of  this  mode.  A 

simple  superposition  of  Figure  4.6.10(c)  on  Figure  4.6.8  leads  to  the  conclusion  that  a 
congregation  of  particles  observed  in  a  still  photograph  does  not  necessarily  correspond  to  a 
concentration  of  vorticity. 

5.  DISCUSSION 

The  velocity  and  length  scales  in  a  plane,  turbulent,  and  supposedly  self-preserving 
wake  are  depend. .nt  on  initial  conditions  and,  therefore,  on  the  shape  and  size  of  the 
obstacle  generating  the  wake.  We  could  not  prove  that  these  scales  will  not  become 
universal  functions  of  fx/G)  at  extremely  large  values  of  (x/9),  but  the  distances  .it  which 
this  may  (or  may  not)  occur  may  be  so  large  as  to  have  no  practical  impact  on  the  problem 
since  the  range  of  x/G  values  in  the  present  study  extended  to  2000  where  uq/U  -  0.03. 
The  dependence  of  the  plane  mixing  layer  on  initial  conditions  was  observed  some  years  ago 
(Champagne  et  al.  1976)  and  traced  to  the  presence  of  large  coherent  structures  (Oster  and 
Wvgnanski  1982'.  The  existence  of  large  coherent  structures  in  a  wake,  hens  ever,  was 
often  confused  with  vortex  shedding,  which  was  so  ably  discussed  by  Karman  (1912)  in  the 
lee  of  a  circular  cylinder  at  low  Reynolds  numbers.  The  large  eddies  proposed  by  Townsend 
(1956)  and  Grant  (1958)  bear  little  resemblence  to  the  structures  observed  presently. 


-32- 


Townsend  suggested  that  the  large  eddies  present  in  fully  turbulent  free  shear  flows  might 
be  generated  by  hydrodynamic  ii. stability  of  the  mean  flow. 

What  is  the  cause  for  the  apparent  dependence  of  the  small-deficit  wake  on  the  shape 
of  the  generator?  The  nature  of  the  flow  in  the  vicinity  of  the  generator,  including  any 
vortices  shed  by  the  generator,  can  provide  the  only  plausible  explanation  for  this 
phenomenon.  The  frequency,  amplitude,  and  predominant  mode  of  the  initial  perturbation 
vary  from  one  geometry  to  another.  For  all  geometries  investigated,  the  predominant  mode 
of  shedding  was  sinuous,  but  the  presence  of  the  varicose  mode  was  also  detected  in  the 
vicinity  of  the  low-solidity  screens  and  circular  cylinder.  The  strongest  sinuous  oscillations 
were  observed  downstream  of  the  thick  symmetrical  airfoil  because  the  initiation  of 
separation  from  one  surface  changed  the  circulation  around  the  airfoil,  moving  the  front 
stagnation  point  toward  the  separated  surface  and  therefore  initiating  a  separation  from 
the  opposite  surface.  The  amplitude  at  the  shedding  frequency  was  three  orders  of 
magnitude  stronger  than  the  background  turbulence.  The  amplitude  of  the  oscillations 
behind  a  circular  cylinder  of  screens  was  approximately  two  orders  of  magnitude  above  the 
background,  while  the  amplitude  of  the  oscillations  downstream  of  the  solid  strip  he!d 
normal  to  the  flow  was  of  the  same  order  of  magnitude  as  the  background. 

It  is  also  suspected  that  the  lower  the  frequency  of  the  shedding,  the  more  persistent 
the  initial  effects  will  be;  in  fact,  the  frequency  of  shedding  downstream  of  the  symmetrical 
airfoil  of  a  given  thickness  depends  on  the  chord  length,  provided  the  flow  separates 
upstream  of  the  trailing  edge.  The  effects  of  frequency,  however,  appear  to  be  less 
significant  than  the  effects  of  amplitude.  A  detailed  investigation  of  the  near  wake  has 
been  undertaken  in  order  to  determine  the  effects  of  initial  condition  more  precisely. 

Cimbala  et  al.  (1981)  observed,  with  the  aid  of  a  smoke  wire,  the  evolution  of  large 
coherent  structures  in  a  wake  of  a  circular  cylinder  up  to  Rej  <  2000  and  in  a  wake  of  two 
screens  at  comparable  Re.  These  structures  became  apparent  some  200  diameters 
downstream  of  the  cylinder  and  had  a  regular  frequency  two  to  three  times  lower  than  the 


Strouhal  frequency.  They  were  unable,  however,  to  corroborate  their  results  with  spectral 
measurements  for  their  high  Reynolds  number  case  and  attributed  it  to  the  high  turbulence 
level  in  their  tunnel.  Some  measurements  of  spectra  at  the  outer  edge  of  the  wake  were 
presently  carried  out  for  two  wake  generators:  (i)  a  circular  cylinder  at  Re0  =  2500 
corresponding  to  Rej  =  5000  and  (ii)  a  screen  having  45  percent  solidity,  also  at  Ree  = 
2500.  In  the  immediate  neighborhood  of  the  cylinder,  the  predominant  spectral  peak  (not 
shown)  corresponded  to  the  shedding  frequency  of  the  cylinder,  i.e.,  at  St^  =  0.206  or  St  a 
=  0.10.  At  x  >  50,  the  predominant  spectral  peak  dropped  quite  abruptly  to  St©  =*  0.03. 
Thereafter,  the  evolution  of  the  spectral  peak  was  rather  slow  and  is  hardly  detectable  on 
the  scale  shown  in  Figure  5.1.  The  spectral  peaks  associated  with  the  screen  tailored  to 
produce  the  same  momentum  thickness  as  the  circular  cylinder  are  similar  to  those 
mentioned  above  at  x  >  250.  The  big  difference  between  the  two  flows  occurs  at  50  <  x  < 
250,  where  the  characteristic  frequency  of  the  spectral  peaks  generated  by  the  screen 
decreases  slowly  with  x.  At  x  >  150,  one  may  detect  the  appearance  of  an  additional  peak 
in  the  spectrum,  which  roughly  corresponds  to  the  spectral  peak  in  the  far  wake  of  the 
circular  cylinder;  this  peak  amplifies  quickly  and  dominates  the  spectrum  at  x  >  ZPO.  It 
seems  that  the  coherent  structures  in  the  near  wake  of  this  particular  screen  retain  some  of 
their  characteristics  up  to  x  =  250,  while  in  the  wake  of  the  circular  cyl  ndor,  tms 
transition  is  accomplished  at  x  =  50.  It  is  believed  that  t lie  shear  'avers  generated  in  the 
wake  of  the  screen  (see  insert  in  Figure  5.1)  undergo  an  instability  process  reminiscent  of 
the  plane  mixing  layer  (Caster  et  a|.  19B4),  generating  eddies  whose  characteristic 
frequency  decreases  in  the  downstream  direction.  Whether  these  eddies  grow  by  entraining 
fluid  from  the  surrounding  stream  or  by  a  process  of  amalgamation  r ..-mains  to  h*  eon. 

Since  Cimbala  et  al.  (1981)  did  not  observe  any  vortex  amalgamation  in  their 
visualization  experiments,  one  would  be  inclined  to  think  that  the  gradual  decrease  in  the 
characteristic  frequency  stemming  from  an  increase  in  the  wavelength  of  these  eddies  is 
caused  by  entrainment.  In  any  event,  once  the  scale  of  these  eddies  becomes  comparable 
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tc  the  width  of  the  screen,  an  interaction  between  two  shear  layers  of  opposing  vorticity 
has  to  occur  before  tPe  wake  will  become  "fully  developed'.  It  is  possible  that  a  chace 
accommodation  ensues,  resulting  in  a  slow  evolution  of  the  typical  large  structures  existmg 
in  the  self-preserving  region.  The  evolution  of  the  fully  developed  wake  is  currently  being 
investigated,  but  the  importance  of  initial  conditions  is  evident  in  Figure  5.1. 

The  mutual  interaction  between  the  large  coherent  structures  and  the  mean  flow  is 
outside  the  scope  of  the  linear  stability  theory.  However,  the  dependence  of  the  mean  flow 
field  on  the  initial  conditions,  and  consequently  on  the  large  coherent  structures,  poses 
precisely  such  a  problem.  The  notion  that  the  free-stream  turbulence  and  the  shape  of  the 
generator  may  have  an  effect  on  the  development  of  a  self-preserving  wake  was  proposed 
by  Symes  and  Fink  (1977).  These  authors  observed  that  the  wake  generated  by  a 
rectangular  cylinder  did  not  evolve  in  the  same  manner  as  the  wake  generated  by  a  circular 
cylinder,  but  the  most  important  observation  was  that  the  evolution  of  the  wake  was 
sensitive  to  grid  turbulence,  whose  integral  scale  was  an  order  of  magnitude  larger  than  the 
scale  of  the  generator.  This  suggested  that  the  externally  imposed  turbulence  interacted 
with  the  flew  far  downstream,  where  the  typical  scales  in  the  wake  and  in  the  free  stream 
became  comparable.  It  also  implied  that  an  instability  mechanism  might  be  responsible  for 
this  result.  The  present  investigation  confirmed  this  notion,  although  nonlinear  terms  should 
be  considered  in  order  to  assess  the  interaction  between  the  mean  flow  and  the  imposed 
oscillation.  Perhaps,  an  integral  approach  similar  to  the  one  used  by  Ko  et  al.  (1970)  might 
predict  such  an  interaction  through  the  Reynolds  stress.  The  concept  of  flow  equilibrium 
and  self-preservation  has  to  be  carefully  reconsidered  in  view  of  the  present  findings,  in 
spite  of  the  fact  that  the  normalized  shape  of  the  mean  velocity  profile  was  not  affected. 
The  dependence  of  the  lateral  distribution  of  u-  on  tin?  njture  of  the  generator  (Figure 
4.3.1)  and  the  relatively  poor  collapse  of  the  dimensionless  data  onto  a  single  G(n)  function 
for  a  given  wake  generator  raise  the  possibility  that  the  flow  is  not  in  equilibrium. 
Although  this  possibility  was  considered  remote  at  the  start  of  this  investigation,  a  piot 
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showing  -he  difference  between  the  maximum  intensity  ( u 2 x  Jncj  :he  intensity  on  the 
centerline  of  the  wake  (u^)CL  normalized  by  Uq  's  shown  in  Figure  5.2  for  the  wakt  of  the 
flat  plate.  If  the  flow  was  in  perfect  equilibrium,  then  (u2)max  -  (u2)q_/u^  should  have 
been  constant  at  all  x.  Although  this  plot  is  very  susceptible  to  experimental  inaccuracies 
and  should  be  treated  with  due  caution,  the  lack  of  constancy  might  have  been  caused  by 
the  modulation  by  the  varicose  mode  of  the  vortices  resuiting  from  the  sinuous  mode  of 
instability.  The  possible  interaction  between  the  two  modes  has  to  be  investigated  in  detail 
by  considering  additional  components  of  the  turbulent  stress  tensor. 

The  assumption  of  parallel  flow  (eq.  3.2.2)  makes  the  eigenfunction  (y)  and  the 
eigenvalues  a  and  B  invariant  with  respect  to  streamwise  distance  from  the  generator. 
Thus,  for  aj  s  0  (i.e.,  spatially  amplified  waves),  only  a  single  mode  containing  the  largest 
|"aj|  need  be  considered.  However,  the  amplification  rates  in  the  small-deficit,  plane 
wake  are  so  small  that  one  cannot  disregard  one  mode  of  instability  in  favor  of  another 

simply  because  its  | -a j  j  is  the  largest.  This  is  so  because  the  long  distances  required  for  an 
unstable  wave  to  amplify  increase  the  relative  significance  of  the  longitudinal  gradients  in 
the  mean  velocity.  By  neglecting  the  varicose  mode  in  favor  of  the  sinuous  one,  Sato  and 
Kuriki  (1961)  had  to  resort  to  nonlinear  effects  in  their  attempt  to  explain  the  cause  for 
the  generation  of  two  rows  of  vortices.  Mattingly  and  Criminale  (1972)  offered  an 
alternate  explanation,  which  is  based  on  the  vorticity  distribution  of  the  sinuous  mode  being 
superposed  on  the  mean  vorticity.  The  generation  of  a  vortex  structure  reminiscent  of  a 
Karman  vortex  street  can  quite  easily  be  attributed  to  the  superposition  of  the  two 
instability  modes,  keeping  in  mind  that  the  most  unstable  frequency  of  the  varicose  mode  is 
only  slightly  higher  than  a  subharmonic  of  the  most  unstable  sinuous  mode.  Sinuous  forcing 
of  the  flow  has  little  effect  on  the  shape  of  the  large  eddies  visualized  by  smoke  until  the 
amplitude  of  the  forcing  becomes  exceedingly  high.  In  this  case,  the  smoke  pattern  is  more 
regular  and  the  eddies  are  located  closer  to  the  centerline  than  in  the  corresponding 


unforced  wake. 
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6.  CONCLUSIONS 


It  was  experimentally  observed  that  the  characteristic  velocity  and  length  scales,  u^ 


and  L  when  suitably  scaled  by  the  momentum  thickness  and  the  free-stream  velocity,  fio 


not  exhibit  universal  behavior  and  do  depend  on  the  initial  conditions  and  therefore  on  the 


geometry  of  the  wake  generator.  The  mean  velocity  profiles  for  each  wake,  when 


normalized  by  their  own  velocity  and  length  scales,  are  self-preserving  and  are  also 


identical  for  all  wake  generators.  The  distributions  of  the  longitudinal  turbulence  intensity 


normalized  in  the  same  manner  are  self-preserving,  but  are  not  identical,  and  depend  on  the 


geometry  of  the  wake  generator. 


Linear  inviscid  stability  theory,  in  which  the  divergence  of  the  mean  flow  was  taken 


into  account,  predicts  quite  well  the  amplification  and  the  transverse  distributions  of 


amplitudes  and  phases  of  externally  imposed  sinuous  waves  in  a  fully  developed  turbulent 


wake  generated  by  a  flat  plate.  It  appears  that  the  large,  coherent,  vortex  structures 


occurring  naturally  in  a  wake  can  be  modeled  by  linear  stability  theory.  Furthermore,  the 


interaction  of  the  two  possible  modes  of  instability  may  be  responsible  for  the  apparent 


Karman  vortex  street-*ype  of  structures  observed  visually  in  the  small-deficit,  turbulent 


wake. 
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FIGURE  CAPTIONS 

1.1  A  sketch  defining  the  nomenclature. 

1.2  Centerline  turbulence  intensity  for  wakes  generated  !jy  circular  cvlir  i  'rs. 

Marasli  (1983),  Re  =  1360,  AR  =  384 
Marasli,  Re  =  3900,  AR  =  123 
A  Townsend  (19? ?),  Re  =  1360,  AR  =  240 
----  Townsend  (1949) 

-  Uberoi  and  Freymuth  (1969),  Re  =  4320,  AR  =  192 

"2^  k  Vamada  et  al.  1980),  Re  =  4000,  AR  =  80 
r”2'  Champ  igne  (1478),  Re  =  19,000,  AR  =  32 
A  Symes  and  Fink  (  1977),  Re  =  6666,  AR  =  150 
A  Symes  and  Fink,  ext.  turb. 

4.1.'  The  shape  of  the  sr  If  similar  nev.  ;  velocity  profile.  Data  for.  0  p  ;rcent 

soll-'ity  '•  c 'e«n  to.  .  *  •  Nio. 

4.2.2  The  variation  of  u0  with  t  for  three  wake  generators. 

Airfoil 

A  70  percent  solidity  screen 
Solid  strip 


4.2.3  The  variation  of  VV^  with  u^/U  for  three  wake  generators. 

OG  Airfoil 

A — A  70  percent  solidity  screen 
(3“^  Solid  strip 

4.2.4  The  dependence  of  9/L  on  Uq/U^.  Symbols  represent  data  from  all  wake 
generators  studied. 

4.3.1  Normalised  turbulence  intensity  distributions  for  three  generators. 

4.4.1  Spectra  of  u  fluctuations  on  centerline  of  flat-plate  wake  at  various  downstream 
positions.  Insert  shows  similarity  of  spectra  when  scaled  by  L  and  U  . 

4.4.2  Spectra  of  induced  u  fluctuations  at  n  =  3  in  the  flat-plate  wake  at  various 

downstream  locations.  Insert  shows  similarity  of  spectra  when  scaled  by  t_0  ar,d 

U  . 

00 

4.4.3  Spatial  amplification  rates  computed  for  several  disturbance  frequencies  in  a 
wake.  Insert  shows  expected  predominant  frequency  (computed)  compared  to 
measured  predominant  frequency  range  as  a  function  of  x. 

4.5.1  The  effect  of  forcing  on  the  measured  u  spectrum  at  x  =  400. 


4.5.2 


The  effect  of  forcing  frequency  on  the  velocity  spectrum  at  a  fixed  location  in 
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4.5.3  Amplitude  distribution  of  u  fluctuations  phase  locked  to  the  external  sinuous 
forcing  signal,  x  =  4C0. 

4.5.4  A  comparison  between  the  measured  and  predicted  u  perturbation  amplitude 
distributions  at  several  x  locations  in  a  sinuously  forced  wake.  A  is  not  A(x)  of 
theory. 

4.5.5  The  variation  of  the  computed  eigenfunction  with  x. 

4.5.6  A  comparison  between  measured  and  predicted  phase  distributions  in  a  wake. 

4.6.1  The  amplitude  distribution  of  the  u  component  f  =  40  Hz  and  x  =  507  in  the 
unforced  wake. 

4.6.2  A  comparison  between  measured  and  predicted  u  amplitude  distributions  at 
several  x  locations  in  an  unforced  wake,  A  is  not  A  x)  of  theory. 

4.6.3  Large  coherent  structures  photographed  in  the*  wake  of  the  flat  plate  with  no 
forcing,  x  =  500  and  Reg  =  600. 

4.6.4  The  spanwise  coherence  of  the  large  eddies  in  the  wake  of  a  circular  cylinder  at 
x  =  600  and  Reg  =  600. 


4.6.5  The  response  of  the  wake  to  various  sinuous  forcing  leveis.  x  =  646  ar A  n  =  0.6. 
O  f  =  50  Hz,  Az/L0  =  10 


0 

f  =  50  Hz,  A  z/Lq  = 

0.4 

□ 

f  =  100 

Hz, 

Az/L0 

=  10 

n 

f  =  100 

Hz, 

Az/L0 

=  0.4 

A 

f  =  150 

Hz, 

Az/Lo 

=  10 

A 

f  =  150 

Hz, 

Az/Lo 

=  0.4 

4.6.6  Coherence  spectra  measured  at  x  =  646,  n  =  0.6,  and  Az:/LQ  =  0.4  in  the  wake  of 
the  flat  plate.  Upper  trace  without  forcing.  Lower  trace  with  forcing  at  50 
Hz. 

4.6.7  (a)  Coherence  spectrum  measured  at  x  =  646,  n=  0.6,  and  Az/L  =  10,  with 
forcing  at  50  Hz. 

(b)  Same  as  (a)  except  n  =  3. 

4.6.8  Streak-lines  computed  assuming  spatial  instability  of  parallel  flow  at  x  =  300: 
(a)  sinuous,  most-amplified  mode;  (b)  varicose,  most-amplified  mode;  (c) 
combined  modes. 

4.6.9  Streak-lines  from  Figure  4.6.8(c)  plotted  on  larger  scale  to  show  staggered 
nature  of  particle  congregation. 

4.6.10  Vorticity  perturbation  contours  computed  for  the  three  cases  of  streak-lines 

shown  in  Figure  4.6.8:  (a)  sinuous,  most-amplified  mode;  (b)  varicose,  most- 

ampiified  mode;  (c)  combined  modes. 
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The  measured  predominant  frequencies  in  the  wake  of: 

(a)  ^  A  circular  cylinder,  Reo  =  2000 

(b)  ilA  The  30  percent  solidity  screen,  Ree  =  2000. 

!JT  5.2  The  variation  of  turbulent  intensity  defect  on  the  centerline  of  the  wake  of  a 

flat  plate. 

V 

,%* 


£ 


'r; 

i 

sr 

*A 

V 


5 


;«% 

V. 


Figure  4.2.2 
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Figure  4.6.4 
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